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1. INTRODUCTION 



Driver safety airbag is under consideration, having an axisymmetric shape when inflated. This 
airbag shape is appropriate for analysis, because a two-dimensional problem may be considered 
instead of a three-dimensional one in this case. This property of the airbag holds if it does not 
interact with other objects. Using an axisymmetric model, we can analyze the airbag shape, pressure 
within, stress state and strength of the airbag shell as a stand-alone mechanical system. 
Though, the purpose of the airbag is to protect a man that sits beside it against impacts and other 
mechanical influences in the case of a car accident. An accident causes parts of the man's body to 
push into the airbag thus increasing its interior pressure and stresses in its shell. The airbag' s shape 
loses its axial symmetry, and the design model becomes three-dimensional and hard to analyze. If 
we knew exact shapes of the man's body which contacts with the shell during an accident and exact 
location of the body, this three-dimensional analysis would have sense. Without knowing the 
shapes, we have to make idealizations somehow, representing the body as a set of orbed objects, 
and define locations of their contact with the airbag. 

One of possible idealizations of the human body's shape is a paraboloid of revolution, and the most 
probable location of its contact with the airbag is the center of the latter. Base of airbag can be 
approximately represented as a plane. In this case, the mechanical model of the system to be 
analyzed becomes more complicated, but the problem remains axisymmetric. 
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2. A SYSTEM OF DIFFERENTIAL EQUATIONS FOR THE AIRBAG ANALYSIS 



2.1. Formulation of the problem 

Determine, by a specific internal pressure within the airbag, the shape of the bag, stresses in the 
shell, the force that pushes the dies off, and the extent of redundancy of the material that causes 
formation of folds. 

The shell of the airbag that is analyzed is made of a soft film of a constant thickness. The film 
has a zero stiffness, both flexural and compression. Under tension, the film behaves as an 
isotropic linearly elastic material. Its relative extensions can be large. 

The shell of the airbag is made in such manner diat it has a shape of a body of revolution without 
folds or stresses in the material when there are no external influences. There can be many shapes 
of this type. One of them is considered to be the initial shape. 

In its operational state, the airbag is filled with gas under some pressure. When deriving 
formulae, we assume the excess internal pressure to be known. 

The airbag interacts with two dies, a flat and a rounded one. The airbag is assumed to lie upon a 
stiff horizontal plane, and the rounded die, a paraboloid of revolution, is located above at some 
distance fi-om the plane (Fig. 1). The axis of symmetry of the airbag is perpendicular to the 
plane. The axis of symmetry of the die coincides with the axis of symmetry of the airbag. Thus, 
the whole system is axially symmetrical. 

When the airbag is inflated, annular folded regions can appear in it because the material does not 
resist the compression. 

2.2. Mathematical formulation 

Let's put the objects of study into a cylindrical system of coordinates r, 9, z so that the axis of 
rotation of the airbag and the rounded die coincide with the axis 

As the airbag is a body of revolution, its shape is completely determined by a fragment of a 
curve that is the meridional cross-section of the airbag's shell (Fig. 1). 
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Fig. 1 .The meridional cross-section of the airbag 



Assume this curve to be defined parametrically using the coordinates r(t), z(t) expressed via the 
parameter /. The parameter t is the length of the curve fragment from the coordinate origin to the 
current point M when the airbag is not deformed. 

Let the fiiU length of the meridional cross-section curve be Tin the initial state, and let the airbag 
have a shape defined by the coordinates ro(0, 2o(0 via the parameter t. 
Let the foundation plane be located at z=0. 

Let the rounded die have the shape of a paraboloid of revolution with the meridional cross- 
section 
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(2.1) +D. 

^ 2R 

The equation shows that the lowest point of the die is at the distance of D from the plane ^=0, 
and the radius of curvature of the die in the vicinity of this point is R. 

To solve the problem, we will use the method of displacements. According to this choice, 
formulae will be derived and written. As all variables are functions of the only parameter t, 
further it is omitted as an argument of the functions. 

2.3. Geometric equations 

Due to the choice of the arc length as the basic parameter, the following relationship holds: 

(2.2) ^|{r:y+uy=u 

Further below, a stroke denotes the differentiation with respect to t. 

The surface of the airbag cannot penetrate to the paraboloid and cannot move down below the 
plane z=0. Therefore, the following conditions must hold: 

the inequality limitation imposed by the paraboloid: 

(2.3) J—^D-z>0. 
2R ^ 

the inequality limitation imposed by the stiff plane: 

(2.4) ij>0 

The conversion of displacements to elastic relative extensions of the airbag material is expressed 
by the formulae: 

(2.5) e,=/-l + Ti,; 

(2.6) E,=--l+i\,; 

where ei is the elastic relative extension of the material in the meridional direction; 
£2 is the elastic relative extension of the material in the circular direction; 
s is the arc length of the meridian of the deformed airbag. . 

(2.7) / = Vff^. 

The dimensionless values r|i and r\2 need additional explanations. These are fractions of the 
material's length consumed by the fold formation. 

T|i is the fraction of the length of meridional fibers absorbed by circular folds; 

r|2 is the fraction of the length of annular fibers absorbed by meridional folds. 

Wavelengths and amplitudes of folds are not the subjects of studies in this report. The folds are 
assumed to be evenly spread across the surface of the aiibag. 

Due to folds being formed at a redundant length of fibers, the following inequalities must hold: 

(2.8) Tii>0; 

(2.9) Ti2>0. 

The application of the formulae ( 2.8 ) and ( 2.9 ) as limitations that express the softness of the 
material is not commonly known. Usually, one uses conditions of nonnegative stresses. The 
usage of these formulae is explained by the method of displacements that will be used. The 
method of displacements must impose limitations upon geometric factors rather than forces or 
stresses. 
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2.4. Physical equations 

Conventional physical equations are under consideration. We assume that the material of the 
airbag is linearly elastic, and its stresses can be expressed through elastic extensions by the 
following formulae: 



where E is the elasticity modulus of the film material; 
V is the Poisson ratio of the film material; 
Gi is the meridional stress; 
02 is the annular stress. 

2.5. Full energy of the system 

We use the variation approach to solution of the problem using the displacement method. This 
approach seems to be reasonable because the problem contains limitations of the unknown 
functions in the form of inequalities. These limitations complicate the derivation of governing 
equations. 

The variation approach formalizes this stage of the problem solution, helps choose a consistent 
system of variables, formulate natural boundary conditions, formulate the conjugate problem if 
necessary [ 2 ], [ 3 ], [ 6 ] Besides, there is no need to look at the pictures to derive the 
equilibrium equations. 

Let's choose the system of fiinctions that would describe the deformed state of the system 
completely. These functions can be r, z, ili, t|2. 

Let's use the variation principle of minimum of the full energy of the system. The full energy, 
including the energy of the load, is expressed by the following ftinctional in this case: 



where h is the thickness of the film; 

p is the excess pressure in the airbag. 

The first term in the brackets represents the elastic energy of the shell's material while the 
second one the energy of gas as an external load. 

The integrand expression should be treated as a fiinction of independent variables r, z, r*, z\'r]u 

Tl2, t. 

2.6. A principle for deriving the equations 

Considering the fact that the problem contains the inequality hmitations ( 2.3 ), ( 2.4 ), ( 2.8 ), ( 
2.9 ) and operates non-smooth (ro) and discontinuous (r'o) functions (in the case of an airbag 
made of two round pieces of film glued together), we will use a pretty general technique known 
as the principle of maximum of Pontryagin [ 4 ], [ 5 ]. This technique was developed to solve 
problems of optimum control, but the Hamiltonian and Lagrangian mechanics can be derived 
from it as special cases. The principle of maximum of Pontryagin operates deftly with 
discontinuous controls and limited regions of variable definition. 

The problem of optimum control consists of finding n state variables Xi(t),..,Xn(t) and r control 
variables ui(t\..,u,(t) in the specified interval of the argument definition ^ ^ ^ that mininiize 
the given functional: 



(2.11) 



(2.10) 



a, =;^-^(e,+veJ; 
l-v 



(2.12) 



n = J (ef + 2ve,e, + )- npr ^ z dt =\f, {r, r, i^^T), ti, ,t)dt . 
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(2.13) =lfo{^y"y\,^U".^ry^)^^''> 

provided the following equations of state hold: 

( 2.14 ) ^' = /(>q,.., \,u^,.,,u^,t) (/=l,..,w). 

The stroke denotes the differentiation with respect to t. 

The following inequality limitations can be imposed upon the control variables: 

(2.15) 2A^,,.;U^)>0 (7=1,..,^). 

Let's introduce n conjugate variables M(0,..,A^/i(0 and build die Hamilton function: 

( 2.16 ) H(>i,.., x,,iVj;,.., iV.,//„..,^,,0 =-/o +Z JV/. . 

/=i 

The principle of maximum of Pontryagin states die following: the optimum control that 
minimizes the functional ( 2.13 ) is achieved at the controls uu..,Ur that satisfy the conditions ( 
2.15 ) and maximize the Hamilton function ( 2.16 ) for each / within the interval to<t< tu 
Also, the following differential dependencies must hold: 

(2.17) H = ~ 0=l,..,n); 

OX. 

(2.18) ^'=1^ (/=l,..,n). 

The relationships ( 2.17 ),( 2.18 ) are canonical equations of Euler for the problem that is to be 
solved. The system of equations ( 2.18 ) is equivalent to the system of state equations ( 2.14 ). 
If the inequality Umitations are imposed upon the variation of the governing variables: 
( 2.19 ) i/^,,, x„,u,,.., u^j)>0 0=lv.,m); 

then the said approach is tme for fragments of the trajectory within a closed region defined by 
the inequalities ( 2.19 ). For fragments of the trajectory located at the boundary of this region, 
one or more limitations ( 2. 19 ) become equalities. 

In this case, the method of Lagrangian coefficients can be used. Let's denote the set of indexes j 
that correspond to the Umitations ( 2.19 ) converted to equalities as J. Let's introduce die 
Lagrangian coefficients \ij(t) (jeJ), and use die function below instead of the function /q: 

(2.20) Fo=/o-ZMy 

This approach will hold also in the case when the expression below is substituted to the equation 

( 2.20 ) instead of all or some 5, 

(2.21) sp=±^fr, 

i=i ox- 

which represent die limitations of the variable variation in die differential form. 
2.7. Formulation of the problem as a problem of optimmn control 

The problem of minimization of the functional ( 2.12 ) at die limitations ( 2.3 ), ( 2.4 ), ( 2.8 ), ( 
2.9 ) can be treated as a problem of optimum control, if we consider the variables r, z to be the 
state variables and r', z\ i\u t|2 to be control variables. At the same time, r\ z! are assumed to be 
independent of r, Let's denote them as Mr, "z- 
The equations of state will be 
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(2.22) r'=ur, 

(2.23) z'=u,. 

The problem includes inequality limitations. Two of them, namely ( 2.3 ), ( 2.4 ), will be 
represented in the form ( 2.21 ), introducing the notation T|3 and TI4 for the obtained functions: 

(2.24) Tl3=-^^.-^ 

( 2.25 ) i\4=Uz. 

This form of the hmitations is more convenient for the subsequent solution of the problem using 
the Runge - Kutta method. 

We'll introduce Lagrangian coefficients into the set of control variables: |i2, II3, ^4, one per 
each limitation, and re-write the function ( 2.20 ) as follows: 

( 2.26 ) Fo +2ve,e2 +^\hnpr^u^-\i,^, -|Li,Ti, . 

Introducing the conjugate variables Nr, Nz, we will obtain the Hamiltonian function: 

( 2.27 )H = -^^(ef +2ve,e, +el)+Kpr'u^ +)IjT1, ^\i,r], +^,ti, +N^u^-hN^u^. 

The Hamiltonian function is assumed to depend upon the following arguments: 
( 2.28 ) H=H(r, z, Nr, N^, J\u ^2, Ur, m„ ^ii, ^2, Ra, m, t)\ 

where t|3 and T|4 can be expressed via the listed arguments by the formulae ( 2.24 ) and ( 2.25 ), 
ei and £2 by the formulae ( 2.5 ), ( 2.6 ). 

According to ( 2.7 ), we will consider s' to be the function of Ur, u{, 
K, h, E, p are constants. 

(2.29) / = 4ul-Vu\\ 
2.8. Derivation of equations 

Using the expression ( 2.27 ), we can obtain systems of differential equations for building any 
fragment of the integral curve, both within the region defined by the inequalities ( 2.3 ), ( 2.4 ), ( 
2.8 ), ( 2.9 ) and belonging to its boundary. Specifying some of the Lagrangian coefficients to be 
permanently equal to zero, we will define a boundary manifold to which the integral curve 
belongs. If some of the Lagrangian coefficients (for example, ^2) is specified to be zero, its 
corresponding limitation (r|2>0) is not obliged to be an equality. To analyze a fragment of the 
curve that belongs to the boundary manifold (t|2=0; T|3=0), we should leave non-zero only |i2 and 
II3 in the Hamiltonian function. The rest of the Lagrangian coefficients and 104) should be 
made zero. If all Lagrangian coefficients were made zero, the Hamiltonian function would 
describe a fragment of the integral curve that lies strictly within the admissible region. 
To shorten writing, let's introduce the following notation: 

(2.30) p,=^; 

(2.31) p,=^; 

( 2.32 ) B=2nrohE; 

(2.33) Ar,=B^^l^; 
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(2.34) jv=B^^^i^. 



1-v^ 

Let's calculate partial derivatives of the Hamiltonian function needed for the subsequent analysis. 
Derivatives with respect to state variables: 

(2.35) -^ = 2nrpu ^-^^u^ ^; 

dr K To 

(2.36) |i = 0. 

Derivatives with respect to conjugate variables: 

(2.38, f=.,. 

Derivatives with respect to control variables: 
(2.39) |^ = N,+fi3P,-i^tiV;; 



(2.40) 



— = N,-\i,+\L, — flV,; 



(2.41) |^ = ^i,-N,; 

(2.42) ^ = ^^-N,; 

Derivatives with respect to Lagrangian coefficients: 

(2.43) |^=n,; 

an 
an, 

dH 

an, 



(2.44) -^ = 11.; 



(2.45) ^ = Tl3; 



(2.46) |^^ = T1,; 

an, 

Using the formulae ( 2.17 ), ( 2.18 ) and the expressions ( 2.35 ), ( 2.36 ), ( 2.37 ), ( 2.38 ), we 
can write the system of differential equations for r, z, Nr, Nz. 

(2.47) r' = ur, 

(2.48) z' = u,; 

(2.49) N' =-2n,pu -^u+^; 

(2.50) N,'=0. 
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The further technique of deriving the equations is as follows. Choose the manifold for which we 
should obtain the differential equations that describe the integral curve. According to 
considerations of Section 2.7, make some Lagrangian coefficients equal to zero. Compose the 
system of equations to determine the control variables and non-zero Lagrangian coefficients on 
the basis of necessary conditions of maximum of the Hamiltonian function, i.e. the conditions of 
its stationarity. To do that, derivatives with respect to zero Lagrangian coefficients must be 
excluded from the system of partial derivative ( 2.39 ) - ( 2.46 ), and the other partial derivatives 
must be assumed equal to zero. 

This system is then resolved with respect to the control variables and non-zero Lagrangian 

coefficients. 



The solutions that have been found are given in Table 1 for six versions of boundary manifolds. 
Table L Formulae for Lagrangian coefficients and control variables 





No folds 


There are folds 


Contact with the upper die 


Vi + P, 


Tli=0 


V^ + Pl 


Tli=0 


p.2=5A 


Tl2=0 


H2=0 


'n2=-A 




K-\A 


1+pf 


K 

«i=- 1 r 

Vi+Pi 


H4=0 


n (k —vA) 


M4=0 


n K 

Vi+Pi 


Free shell 




ili=0 




Tli=0 




Tl2=0 


\L2=0 


Tl2=-A 


^3=0 


_nXk-va) 


^3=0 




M4=0 


_N,{K-vA) 


M4=0 


»2 = 1 


Contact with the lower 
die 






\ll=Nr 


tli=0 


\i2=BA 


Tl2=0 


\l2=0 


T|2=-A 


H3=0 


wi=^-vA 


^3=0 


ui=K 


\U = -No 


1*2=0 


M4=-A^o 


«2=0 



This table uses the following notation (also, further below): 

(2.51) No=N,+npr^; 

(2.52) N^=piNr-No; 

(2.53) ^:=1+^; 

B 
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(2.54) A = p,-1+^. 



To determine K and A, the value of )ii should be taken from its own boundary manifold. 
Table 2. Systems of differential equations for state and conjugate variables 





No folds 




niere are folds 


Contact with the upper die 




Transition: A>OT; A<Oi 


K 

^ " 1 T 

Vi+pi 




^- Pl^ 

Vl+pf 


BA N. ^K-vA 
N= + 2nrpp, ^ — , 






t (i+pf)RJVi+pf 




n;=o 




Transition: A^a^OT; NA<Oi 




Transition: N/;>Ot; NakOI 


Free shell 


, N^K-vA) 
^ ^ 1—. 7 


Transition: A>OT; A<Oi 




,_iVo(K-vA) 




^^_BA 2nrpN,{K-vA) 


= 1 


ISt =0 


AT =0 




Transition: A/i^T; A/b<Oi 




Transition: A^o^t; A/o<Oi 


Contact with the lower die 


r'=K-vA 


Transition: A>Ot; A<Oi 




z'=0 






iv;=o 


n;=o 


AT =0 



Substituting the relationships given in Table 1 to the formulae ( 2.47 ) - ( 2.50 ), we will obtain 
the systems of differential equations of interest. They are given in Table 2. 



This table lists also conditions of transition from one boundary manifold to another. Three 
determinant fimctions are used to do that: A^o, A^a and A. Alteration of sign of the A^o function 
shows the separation of the airbag and the flat die. Alteration of sign of the A^a function shows 
the separation of the airbag and the rounded die. Alteration of sign of the A function occurs at the 
boundary of the smooth and folded regions of the airbag. These functions are the determinant 
ones for Lagrangian coefficients. 
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We suggest to use the principle of extension by parameter. If the boundary manifold has been 
determined correcdy in the initial point of integration, then the transition from one boundary 
manifold to another, according to altered signs of the determinant functions, will have the correct 
result. 

Note that all the relationships have been obtained without any simplifications. Both 
displacements and relative extensions can be arbitrarily large. 

2.9. Boundary conditions 

The boundary conditions for the problem to be solved can be formulated in a pretty simple form: 
( 2.55 ) r(0)=0; z(0)=0; 

(2.56) r(T)=0',z(T)=D. 

Though, these boundary conditions are inconvenient for use with the methods of Euler or 
Runge - Kutta. It would be better to know values of all four unknown functions at the same point 
of the integral curve. 

The second inconvenience is that the process of integration can become unstable when the 
integral curve goes from "soft" folded regions to "stiff" regions without folds, or to regions of 
contact with the die. 

To eliminate this second inconvenience, we suggest the technique of integration "from edges to 
the middle". According to this technique, two integral curves are searched for using the Euler or 
Runge - Kutta method. One of them is built for the range of t values [0..r/2]. It begins in the 
initial point and ends in the middle point of the integration interval. The second curve is built for 
the range of t values [T.,T/2], It starts in the initial point and ends, just the same as the second 
curve, in the middle point of the integration interval. All four unknown functions r, 2, Nr, Nz must 
be consistent in the point r=r/2. 

To begin the integration, we have to know also values of conjugate variables A^;^ and A^^, besides 
knowing the values ( 2.55 ), ( 2.56 ). Four values of the conjugate variables at the ends must be 
selected in such manner that the four unknown functions r, z, A^^, Nz be consistent at t=T/2, The 
function Nz is constant on the whole interval of integration, and this facilitates the task to some 
extent. Its values at the ends of the integration interval are equal: Nz (0)=Nz{T), 

There is another difficulty to be mentioned. In the two edge points of the integral curve, at r=0, 
z=0 and r=0, z=D, some auxiliary functions become infinite. Therefore, the integration is better 
started at some gap from these points. To be definite, let this gap be At. 

Assuming that the material of the shell is evenly tensioned in the vicinity of the point r=0, we 
can state for this vicinity: 

(2.57) X=l+-^; 

1-v 

(2.58) A = — ; 

1-v 

(2.59) K-vA=l-fe. 

Then, using the systems of differential equations included in Table 2, we will obtain initial 
conditions for the point defined by the parameter t=T-At: 

(2.60) r(r-A/)=A/(l + e); 

(2.61) ^(T-A/) = D+M1; 

(2.62) Ar^(r-A/)=-A/[^..i^]; 
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(2.63) N^(r-A/) = N^. 

In the same way, we will obtain the following for the point at t=At: 

(2.64) r(A/)=A/(l + e); 

(2.65) ij(A/)=0; 

(2.66) N,(A/)= ; 

(2.67) iV(A0=N^. 

In the point of contact with the upper die, e is assumed to have its particular value. In the point of 
contact with the lower die, another value of e is assumed. Though, values of A^^ in those points 
must be equal. 

3. LAWS OF DEFORMATION OF AXISYMMETRIC AIRBAGS 
3.1. Assumptions 

The design of the system of interest is shown on Fig. 1. 
The notation: 

Rq is the radius of the aiibag in the non-deformed flat state; 
h is the thickness of the film that the shell is made of; 
E is the film's elasticity modulus; 

V is the Poisson ratio; 

r, 9, z are the coordinates, namely: 
r is the radial coordinate; 
0 is the angular coordinate; 

z is the Unear coordinate counted along the rotation axis; 

t is the length of the line that connects the coordinate origin and the current point M of 
the meridional cross-section of the airbag in the non-deformed state; 

ro is the radial coordinate for the non-deformed airbag; 

R is the curve radius of the moving body; 

D is the distance between the moving body and the support plane; 

Q is the force of pressure of the moving body upon the airbag; 

p is the overpressure within the airbag; 

V is the volume of the airbag in the deformed state; 

01 is the meridional stress; 

02 is the circular stress; 

T|2 is the fraction of the length of circular fibers absorbed in meridional folds. 

As the airbag is made of two flat round sheets of a material of the constant thickness of the 
radius Rqj the following dependency holds and is used in computations: 

(31) 

ro=2Ro-/ (/>Ko) 
3«2. The basic design of the airbag 

Numerical computations have been done using the units of kilogram-force (kg-f) and centimeter. 
The following basic airbag design parameters were agreed upon: 
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( 3.2 ) /?o=35 cm; h=Om cm; £=35000 kg-£^cm^; v=0.39; R=40 cm;p=0.2 kg-^cml 

Fixed-parameter calculations were done exactly with the basic design. When doing the serial 
analysis within a certain range of variation of tfie parameters, the basic design was the central 
point. 

3.3. The dependency between the stress and strain state of the airbag and the material 
rigidity 

A series of calculations was performed for the airbag that did not touch any bodies, with the 
following source data: 

( 33 ) /?o=35 cm; /i=0.01 cm; v=0.39;/7=0.2 kg-f^cml 

The elasticity modulus of the material varied in the range of £=[63 6.. 3500000] kg-^cm^. The 
lower boundary of this range is close to the critical value of 593.5 kg-fi'cm^ at which the 
displacements rise dramatically even at a slight decrease of the elasticity modulus. The upper 
boundary is 100 times greater than the basic value. 

Results of the analysis are represented by charts. The shapes of the airbags depending upon the 
material's elasticity modulus are shown on Fig. 2 and Fig. 3. 

Fig. 2 shows meridional cross-sections of the airbags made of soft films with a low elasticity 
modulus similar to rubber. At these elasticity modulus values, the airbags inflate immensely, 
acquiring nearly spherical shapes. Folded regions disappear. The surface of the airbag becomes 
smooth. Relative extensions of the material are large. 

70 
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E=3500 

Fig. 2. The airbag shape in the case of a pliable film. /?o=35; /i=0.01; v=0.39; /?=0.2 
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E=11068 

E=35000 

£=350000 

E=3500000 



Fig. 3. The airbag shape in the case of a rigid film. /?o=35; /i=0.01; v=0.39; p=02 

Fig. 3 shows meridional cross-sections of the airbags made of rigid films. The elasticity modulus 
of the shell film is 10000 kg-f^cm^ to an almost infinite value. In this series of computations, the 
airbag shape depends httle upon the material's elasticity modulus. The airbag looks like a 
flattened ellipsoid of revolution with the aspect ratio 0.6:1. There are intensive folds formed in 
the equatorial region of the airbag. 
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Fig. 4. Regions and intensities of fold formation. /?o=35; /i=0.01; v=0.39; /7=0.2 

Fig. 4 shows regions where folds are formed and diagrams of the fraction of the circular fiber 
length of the airbag absorbed by the meridional folds. The argument is the t parameter, the length 
of the line that goes from the coordinate origin to the current point of the meridional cross- 
section of the airbag in the non-deformed state (see Fig. 1). The value of the parameter r=0 cm 
corresponds to the center of the round sheet that lies on the support plane, r=35 cm is the glued 
seam of the round sheets, t-lQ era is the center of the sheet that faces the man. 
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Fig. 5. Meridional stresses for the pliable film. /?o=35; /z=0.01; v=0.39; p=02 

As you can see from the diagrams, the more rigid the shell material is, the wider area is occupied 
by the folds, and the more fraction of the circular fiber length is absorbed by the folds. At the 
elasticity modulus of £=636 kg-f^cm^ and £=1107 kg-fi^cm^ the folds do not appear at all, and 
the surface of the airbag becomes smooth. 

The stress state of the airbag's shell can be described by meridional and circular stresses. Fig. 5 
shows diagrams of meridional stresses in a pliable film, while Fig. 6 shows those in a rigid film. 




0 5 10 15 20 25 30 35 40 45 50 55 60 65 70 
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E=l 106800 

E=3500000 

Fig. 6. Meridional stresses in a rigid film. Rq=35] /i=0.01; v=0.39;p=0.2 

As the diagrams show, the maximum stresses are found in the central areas of the round sheets. 
The minimum stresses appear in the seams. 
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Fig. 7. Circular stresses in a pliable film. /?o=35; /i=0.01; v=0.39; p=02 

The distribution of circular stresses is shown on Fig. 7 and Fig. 8. The circular stresses in the 
central areas of the shell are equal to the meridional stresses. At a distance from the central area, 
both circular and meridional stresses decrease. The circular stress decreases faster. 
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Fig. 8. Circular stresses in a rigid film. Ro=35\ A=0.01; v=0.39; p=02 
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E/35000 



Fig. 9. Dependence of the airbag thickness D on the elasticity modulus. 
/?o=35; /i=0.01; v=0.39; p=0.2 

Now let's discuss the diagrams of dependence between the stress and strain parameters of the 
airbag and the film's tension rigidity, shown on Fig. 9 - Fig. 12. These diagrams include the 
relative elasticity modulus of the material along the abscissa axis. The unit of measurement is the 
elasticity modulus of the basic airbag design ^=35000 kg-^m^. The logarithmic scale is used 
here due to a wide range of the material's rigidity variation (almost 10,000 times). 
All diagrams include the critical left point that corresponds to the minimum elasticity modulus of 
£=593.5 kg-f^cm^. At this point, the stable mode of equilibrium of the airbag changes to the 
instable mode. From the physical viewpoint, it means that the further decrease of the elasticity 
modulus at a specific internal pressure will result in the avalanche-like rise of deformations and 
the destruction of the shell. 

Fig. 9 shows the diagram of the dependence of the airbag thickness D on the elasticity modulus. 
At large values of the elasticity modulus, the airbag thickness varies litde and is about 32 to 
33 cm. As the film rigidity falls down, the airbag inflates to the thickness of 90 cm, acquiring the 
spherical shape, and then it breaks whatever the film strength might be. 




Fig. 10. Dependence between the airbag volume Vand the elasticity modulus. 
/?o=35; /i=0.01; v=0.39; /i=0.2 

The volume of the airbag that is studied (Fig. 10) made of a rigid film is 52.3 liters. The critical 
volume of the airbag is 488 liters. 
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Fold formation in the place where the round sheets are glued (Fig. 11) is intensive and absorbs 
about 24% of the material of the rigid film. As the rigidity goes down, folds become less clear 
and intensive, and they disappear at the elasticity modulus of 1 150 kg-^cm^. 

0.25 



0.2 



0.15 
0.1 



0.05 



%.0l 0.1 1 10 100 

E/35000 

Fig. 11. Dependence of the fold formation intensity on the elasticity modulus. 
/?o=35; /i=0.01; v=0.39; p=02 

Fig. 12 shows the dependence of stresses in the bag material vs. the elasticity modulus. The solid 
line shows meridional and circular stresses (equal to each other) in the central area of the bag. 
These stresses show minimum values at the elasticity modulus of about 2000 kg-:^cm^. The 
minimum value is 456 kg-f^cm^. As the film rigidity falls, the stresses rise to the critical value of 
1300 kg-fi^cm^. As the rigidity rises, these stresses also rise infinitely. 
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Fig. 12. Dependence of the stresses in the bag material on the elasticity modulus. 

/?o=35; /i=0.01; v=0.39; p=0.2 

The dotted line shows the meridional stresses in the glue seam of the film sheets. These stresses 
are not caused to increase by an increased film rigidity, and they tend to achieve 203 kg-^cm^. 
Their critical value is 744.5 kg-f^cm^ at the film rigidity of 593.5 kg-^cm^. 

3.4. Deformation of the airbag by the human body 

This section studies the interaction between the airbag and the human body for the basic design 
of the airbag. The human body is simulated by a paraboloid of revolution (Fig. 1) with the 
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curvature radius in the vertex R=40 cm. Depending on the force caused by the man's pressure 
upon the bag, the airbag's shape, its stress and strain state, its internal volume and the depth of 
the man's body indentation are determined. 

The shape of the airbag is shown on Fig. 13 and Fig. 14. The analysis was done with fixed values 
of the interaction force between the man and the airbag Q. The step of variation of this force is 
50 kg-f, except for the last value of 23.5 kg-f that corresponds to a completely flat airbag. 

Fold formation in all phases of flattening can be observed in the equatorial region of the airbag. 
At an intensive flattening, the folds appear even where the man's body touches the bag. The 
regions and intensities of fold formation are shown as diagrams on Fig. 15 and Fig. 16. 

Just as in the case of a free airbag, the stresses show their minimum in the equatorial region, and 
they rise when approaching the poles. The maximum values of the stresses are found in the 
center of the airbag, exactly where it touches the man. The distribution of meridional stresses is 
shown on Fig. 17 and Fig. 18, and the distribution of circular stresses is shown on Fig. 19 and 
Fig. 20. Also, Fig. 21 shows diagrams of stresses in the central part of the bag vs. degree of 
flattening of the airbag. 
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Q=100 

Q=150 

Q=200 

Fig. 13. Shapes of the airbag at a slight flattening 
/?o=35; /?=40; /i=0.01; £=35000; v=0.39;p=0.2 
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Fig. 14, Shapes of the airbag at an intensive flattening 
/?o=35; /?=40; /i=0.01; £=35000; v=0.39;p=0.2 
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Fig. 15. Regions and intensities of fold formation at a small flattening 
/?o=35; /?=40; /i=0.01; £^35000; v=0.39; p=0.2 
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Fig. 16. Regions and intensities of fold formation at a large flattening 
/?o=35; /?=40; /i=0.01; £=35000; v=0.39;p=0.2 
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Fig. 17. Meridional stresses in the airbag at a small flattening. 
/Jo=35; /f=40; /i=0.01; E=35000; v=0.39;p=0.2 
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Fig. 18. Meridional stresses in the airbag at a large flattening, 
/?o=35; /?=40; /j=0.01; £=35000; v=0.39;/7=0.2 
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Fig. 19. Circular stresses in the airbag at a small flattening. 
/?(P35; /f=40; /i=0.01; £^35000; v=0.39;/?=0.2 
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Fig. 20. Circular stresses in the airbag at a large flattening. 
/?o=35; /?=40; /i=0.01; £=35000; v=0.39; p=02 
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.21. Stresses in the central areas of the bag vs. the degree of the bag flattening. 
/?o=35; R=40] h=Om; £=35000; v=0.39;/?=0.2 
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Fig. 23. The volume of the aiibag V vs. the degree of the bag flattening. 
/?o=35; R=40; /i=0.01; £=35000; v=0.39; p=0.2 

It is much of interest to know the dependence of the airbag's volume and its reaction force on its 
degree of flattening. These data are necessary to solve the dynamical problem of softening the 
man-car impact during the car collision with a barrier. The mentioned dependences are shown on 
Fig. 22 and Fig. 23. 

3.5. Dependence of the stress and strain state of the airbag on its internal pressure 

The assumption that the strains depend Uttle on the intemal pressure while the stresses are 
proportional to that seems pretty reasonable. To validate this assumption, a series of problems 
have been solved at the following fixed initial data: /?o=35; /f=40; /i=0.01; £=35000; v=0.39; 
D=17.5. The range of the intemal pressure variation is /7=0.02..2 kg-G'cm^, that is, the topmost 
boundary is one hundred times the lowest pressure and the basic value is 0.2 kg-f^cm^. 

Results of this study are shown on Fig. 24 to Fig. 26. They confirm the mentioned assumption 
only partially. 

The abscissa axes of all diagrams shows the relative variation of intemal pressure within the 
airbag p/po where po=02 kg-fi^cm^, i.e. the pressure assumed for the basic design. The ordinate 
axis shows values of various functions normalized for their values in the basic design. 
The study shows (see Fig. 24) that, indeed, the dimensions of the airbag change Utde as the 
pressure varies. If we exclude too great a value of p=2 kg-f^cm^ then the radial coordinate will 
vary within the Umits of -1% +2%. The volume of the airbag experiences slightly larger 
variations, within the Umits of -2% .. +5%. 
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Fig. 24. Relative variation of the airbag volume V and the radial coordinate r(35) at the glue 
seam vs. relative pressure within the airbag p 

As for force/stress factors, their relative changes look like straight lines on diagrams in the 
logarithmic scale. It means that the dependence like that can be described by the following 
formula: 




(3.4) 



where F is the value of the unknown function at the pressure p\ 
Fq is the value of the unknown function at the pressure poJ 
a is an empirical value. 
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Fig. 25. Relative reaction force of the airbag Q vs. relative pressure within the airbag p 
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Fig. 26. Relative stresses vs. relative pressure in the airbag p 
4. CONCLUSIONS 

• The efficiency of the described technique for analyzing the airbag has been confirmed by a 
lot of practical solutions obtained. 

• Results of the analysis show that the film of the thickness /i=0.01 cm experiences 
unallowably large stresses. It is necessary to increase the thickness of the film or (and) its 
strength. 

• Dependencies obtained in this report enable one to adjust properties of the airbag and choose 
a material for its fabrication. 

• The dependences of the reaction force Q of the airbag and its volume on its degree of 
flattening D have been obtained. These dependences enable one to pose and solve a 
dynamical problem of interaction between a man and the airbag after the car impacts an 
obstacle. The solution of this dynamical problem will enable us to trace the following values 
in time: internal pressure within the airbag, the force of interaction between the man and the 
airbag, the velocity acquired by the man when pushed off by the airbag, and other physical 
parameters. 
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